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■ Abstract. The physical reasons why the Drudc dielectric function is not compatible 

with the Lifshitz formula, as opposed to the generalized plasma-like permittivity, are 
presented. Essentially, the problem is connected with the finite size of metal plates. It is 
shown that the Lifshitz theory combined with the generalized plasma-like permittivity 
is thermodynamically consistent. 
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CN ■ 1. Introduction 

q 

In the last few years the Casimir effect [Tj received common recognition as one of the 
\ most important subjects of interdisciplinary interest. The Casimir force is of the same 
nature as other one- loop vacuum effects of quantum electrodynamics [2J. It arises due 
to the alteration of the spectrum of electromagnetic zero-point oscillations by material 
boundaries. Early stages of modern experiments and related theory are reflected in [3J. 
Recent trends go towards complex experimental and theoretical studies of the Casimir 
effect, including the applications to nanotechnology. For this purpose many classical 
theoretical results on the subject obtained within the framework of quantum field theory 
(see, e.g., monographs [UG3E], reviews [3], [HE] proceedings [9] and more recent papers 
[10] [XT] ) should be adapted to the case of real material bodies. Realistic material 
properties are important also for applications of the Casimir effect in nanotechnology 

(Tansies]. 

The basic theory of the van der Waals and Casimir forces between dielectric 
materials was proposed by Lifshitz [15l . However, the application of this theory to 
Drude metals and semiconductors with sufficiently low charge carrier density met serious 
problems. Namely, at first, it was shown [161 [TT] that for Drude metals with perfect 
crystal lattices the Lifshitz theory violates the Nernst heat theorem. This is connected 
with the fact that the reflection coefficient for the transverse electric mode of the 
electromagnetic field at zero frequency is equal to zero if the dielectric permittivity 
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behaves as when the frequency uo vanishes. For metals with impurities the Nernst 
heat theorem is formally preserved [18J. This inclined the proponents of the Drude model 
to believe that it is applicable together with the Lifshitz theory (different arguments on 
this problem can be found in [191 EH EH [22]). However, precision measurements of the 
Casimir force [23], EH E5] [26] excluded the Drude model at a 99.9% confidence level. 

Another problem arises when the Lifshitz theory is applied to dielectrics or 
semiconductors with not too high density of free charge carriers. In this case the Nernst 
heat theorem is violated if the conductivity at zero frequency is taken into account 
[27] [28] [29l [30] . Here, it is the discontinuity of the transverse magnetic mode at zero 
frequency, which is responsible for that violation. Furthermore, it was demonstrated 
[31] [32] that the inclusion of the conductivity at zero frequency into the Lifshitz theory 
leads to a contradiction with experiment. Until recently there exists not any theoretical 
approach to the thermal Casimir force which would be in agreement with both short 
separation experiments [33] [33] and longer separation experiments [23] E3] [25] [26] . The 
approach using the usual, nondissipative, plasma model was shown to be in agreement 
with longer separation experiment [23] EU [25] [26], but to be in contradiction with the 
experiment performed at short separations [33], [33]. The impedance approach [35J was 
also found in agreement with longer separation experiments [23] EU E5] [26], but it is 
simply not applicable at shorter separations characteristic for the experiment of [33] [33]. 
Because of this, the Lifshitz theory at zero temperature by necessity was used for the 
comparison between the short separation experiment [33] [33] and theory, even though 
that experiment was performed at a room temperature of 300 K. 

Recently, a new theoretical approach to the thermal Casimir force between real 
metals has been proposed [36J using the generalized plasma-like dielectric permittivity. 
The latter includes dissipation processes due to the interband transitions of core 
electrons but disregards dissipation due to scattering processes of free electrons. As was 
shown in [36], the Lifshitz formula combined with the generalized plasma-like dielectric 
permittivity is consistent with both short and long separation experiments. It also 
exactly satisfies the Kramers -Kronig relations. However, the question why one should 
include one type of dissipation (interband transitions of core electrons) to fit theory 
to experiment while disregarding another one (scattering processes of free electrons) 
remained unresolved. 

In this paper we present and discuss the physical explanation why the Drude 
dielectric function cannot be used to describe the thermal Casimir force between metal 
plates of finite area. The idea of that explanation was briefly published first by Parsegian 
[37] . but did not attract the attention it is deserving. As we show below, the Drude 
dielectric function is not compatible with the zero-frequency term of the Lifshitz formula 
if the area of plates is finite. We also perform a rigorous analytical proof of the fact that 
the Casimir entropy calculated using the generalized plasma-like dielectric permittivity 
satisfies the Nernst heat theorem. 

The paper is organized as follows. In Section 2 we explain why the Drude dielectric 
function is incompatible with the Lifshitz formula in the case of two parallel metallic 
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plates of finite area. Section 3 contains an asymptotic derivation of the analytic 
expression for the Casimir entropy in the limit of low temperatures. Here we present a 
proof for the validity of the Nernst heat theorem in the Lifshitz theory combined with 
the generalized plasma-like dielectric permittivity. Section 4 contains our conclusions 
and a discussion. 



2. Why the Drude dielectric function is not compatible with the Lifshitz 
formula for metallic plates of finite area 

In the framework of the Lifshitz theory the free energy of the fluctuating electromagnetic 
field between two electrically neutral plane parallel plates of thickness d at temperature 
T in thermal equilibrium is given by [3], [15] 

k R T ~ / 1 



F(a, T) = 22- 53 (l - ^oi) r k ± dk ± 

2lX JZq V 2 J Jo 
x {in [l - 4 M (6, fe ± )e- 2a «] + In [l - r^fe, ^)e" 2 ^] } . (1) 



Here a is the separation distance between the plates, ks is the Boltzmann constant, 
= 2irkBTl/h are the Matsubara frequencies defined for any I = 0, 1, 2, . . . , and 
k±_ = \k±\ is the magnitude of the wave vector projection onto the plane of the plates. 
The reflection coefficients for the two independent polarizations of the electromagnetic 
field (transverse magnetic, TM, and transverse electric, TE) are expressed [38] in terms 
of the frequency-dependent dielectric permittivity, e(u), along the imaginary frequency 
axis: 

eh? - kf 



ttm(€i, k 



± 



e iQi + kf + 2qihei coth(M) 



kf - qf 

rTE& ' k±) = qf + kf+^kcoMkdV (2) 



where 



^V^i + f' k^^jkl + e^, e, = e(iei)- (3) 

Equation ([T]) was originally derived [15J for dielectric plates of infinite area. 
However, it is commonly used for plates of finite area S under the condition a <C V^- 
If this condition is satisfied, corrections to Eq. (CQ) due to the finiteness of plate area are 
shown to be negligibly small [3j [39] for both dielectric and ideal metal plates. Below 
we show that this is not the case for metal plates described by the Drude dielectric 
function where the presence of a real current of conduction electrons leads to a crucially 
new physical situation. 

Papers [HI EH [22j HQ] describe metallic plates by using the dielectric permittivity 
of the Drude model, 

eD{«>) = 1 - -4-, (4) 
uj(<jJ + 17) 
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where uj p is the plasma frequency and 7 is the relaxation parameter. As is correctly 

stated by Parsegian (see [37], p. 254), "this is valid only in the case of an effectively 

infinite medium where no walls limit the flow of charges." To gain a better understanding 

of this statement, we derive Eq. starting from Maxwell equations in an unbounded 

nonmagnetic metallic medium 

„ IdE 4vr „ „ 
rotB = -— + — a E, dwB = 0, 

c at c 

1 dB 

rotS = — , divE = 0. (5) 

c at 

Here, the electric current density j = a Q E is induced in a metal under the influence of 
external sources, and a is the conductivity at zero frequency. Physically the demand 
that the medium is unbounded means that it should be much larger than the extension 
of the wave fronts of electromagnetic waves coming from external sources (i.e., of zero- 
point oscillations and thermal photons). 

Solutions of equations §J5§ can be found in the form of monochromatic waves, 



E = Re 



E (r)e-^ , B = Re B (r) e -^ , (6) 



where E (r) and B (r) satisfy equations 

AE (r) + k 2 E (r) = 0, AB (r) + k 2 B (r) = 0, (7) 
following from (jSJ) with 

2 uj 2 A7ia uj _ e n (uj)uj 2 

Here the dielectric permittivity of the normal skin effect, s n (u), is introduced 

e n ( w =l + i °. (9) 

UJ 

This equation is applicable at not too high frequencies (the region of the normal skin 
effect) where the relation j = a E is valid. The Drude model extends the applicability 
of ([H]) to higher frequencies, up to the plasma frequency, by making the following 
replacement in equation Q: 



(To 1 + - 

a -> a(u) = } ■ J ' . (10) 
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Substituting fflUj) in (jUJ) and taking into account that a = uj 2 / (Attj) [41J we recover 
the dielectric permittivity of the Drude model (T4|). At sufficiently high frequencies 
7 uj < uj p (the region of infrared optics) one can neglect unity as compared to cj/7 
and uj 2 /'-/ 2 in (flOl) . Then (T4j) and (TTTJT) lead to the so-called free electron plasma model 

e p = l--^, a{uj) = . 11 

U3 l UJ 

Thus, the plasma model is characterized by pure imaginary conductivity. In the opposite 
limit w<7 the unity in both numerator and denominator of (TTOT) dominate over cj/7 
and uj 2 /7 2 leading to <j(u) = cr . This converts the dielectric permittivity of the Drude 
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model dlj) in the dielectric permittivity of the normal skin effect characterized by 
real conductivity of conduction electrons a . 

The total current in the framework of the Drude model (jlj) is given by 

iuj 



3u>t( r ,t) = Re 



-—e D {u)E (r)e 



-ihjt 



(12) 
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The first term on the right-hand side of this equation has the meaning of the 
displacement current, whereas the second term, in accordance to (J6j), is proportional 
to the physical electric field E = E(r,t) and describes the real current of conduction 
electrons. Under the condition 7 < w < u p , i.e., in the region of infrared optics, the 
first term dominates. This is the displacement current of the plasma model with a pure 
imaginary conductivity ( Tl2l . Under the opposite condition uj <C 7, i.e., in the region 
of the normal skin effect, the second term on the right-hand side of ffl2|) . i.e., the real 
physical current of conduction electrons dominates. 

After the above discussion on the derivation of the Drude model we now return 
to the role played by the finite size of the plates. Let us consider plane waves of zero- 
point oscillations and thermal photons in between the plates incident on their interior 
boundary surfaces. It is common knowledge (see, e.g., [121 H3]) that charge carriers in 
a conductor are moving in response to electromagnetic oscillations. For plates of finite 
area the application condition for the derivation of the Drude model is formally violated 
because the extension of the oscillation wave fronts is much larger than the size of any 
conceivable plates. However, if the frequencies of oscillations are high enough (recall that 
at room temperature the first Matsubara frequency is equal to ^ « 2.47 x 10 14 rad/s, 
and all others with I > 1 are, respectively, higher) there is no accumulation of charges 
on the side boundary surfaces of finite metal plates. First, at so high frequencies the real 
current of conduction electrons [given by the second term on the right-hand side of ([12"]) ] 
is small in comparison with the displacement current. Second, a high-frequency electric 
field quickly changes its direction. As a result, electric charges which are accumulated 
on the sides of a plate change their sign many times during any reasonable time of force 
measurement. This leads to practically zero mean surface charge. Thus, equations (jlj), 
(J2J), jE]) and (TT2]) remain macroscopically valid. 

The situation changes drastically when the contribution from the zero Matsubara 
frequency £0 = is considered. The plane wave of zero frequency should be understood 
as limit of plane waves with some low frequencies £ in the case that £ — > 0. As was 
mentioned above, the extension of a wave front far exceeds the size of the plates. If 
it is remembered that the period of the wave of vanishing frequency goes to infinity, 
the Casimir plates are found in practically constant electric and magnetic fields. 
As is described in textbooks on classical electrodynamics (see, e.g., (121 S3]), in a 
quasistatic case the propagation direction of a plane wave inside a metal is approximately 
perpendicular to its surface independently of the angle of incidence. Thus, a short-lived 
current which arises under the influence of a constant electric field in the plane of 
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plates immediately gives rise to some nonzero surface charge densities p of opposite 
signs accumulated on opposite sides of the plates. The electric field generated by these 
charges precisely compensates the electric field of zero frequency inside a metal. As a 
result, the electric field inside the metal is exactly equal to zero [121 [43]. As to the space 
between the plates, the resulting field there is the superposition of an approximately 
constant field of external sources and of the field generated by the charge distribution 
on the plate sides. 

From what has been said, it appears that Maxwell equations (jSJ) and all 
consequences obtained from them are not applicable in the case of plane waves of 
zero frequency. In that case for finite plates not only a nonzero induced current must 
be taken into account but also a nonzero induced charge density generated by this 
current which, however, is omitted in (JSJ). As was noted by Parsegian (see [37J, p. 254), 
"conductors must be considered case by case corresponding to the limitations imposed 
by boundary surfaces." Here we demonstrate that these limitations arise from the 
substitution of the Drude dielectric function (which is obtained for unbounded medium) 
in the zero-frequency term of the Lifshitz formula. Such substitution is in contradiction 
with electrodynamics, because, as was shown above, the Drude model admits a nonzero 
current of conduction electrons, whereas electrodynamics assertains that the current of 
conduction electrons inside a finite metal plate placed in a plane wave of zero frequency 
must be equal to zero. 

We emphasize that the time interval during which charges on the sides of finite 
plates are accumulated and the total electric field in a metal turns into zero is extremely 
short. To make sure that this is the case we describe the dynamic process of charge 
accumulation on the sides of the plates by a simple model: 



where E is the constant electric field along the plates and E p (t) is the field produced by 
the surface charge density till the moment t. For the order of magnitude estimation it 
is sufficient to represent E p (t) as the field in a plane capacitor: E p (t) = 4irp(t). Then 
we arrive at the solution 



As an example, for Au it holds ArcaQ = 3.5 x 10 18 s -1 and, thus, even after a 
very short time lapse of t = 10~ 18 s, p(t) practically achieves the maximum value 
p(oo) = E/(4n). Then from equation (TT3"]) it follows that the total field inside a metal 
vanishes, E tot (oo) = 0, as it should be. 

Thus, the substitution of the Drude dielectric function in the zero-frequency term of 
the Lifshitz formula is self-contradictory. As was recalled in Introduction, for the Drude 
model it holds ?"te(0, k±) = 0. This means that the TE field of zero frequency completely 
penetrates into a metal. For metal plates of finite size this inavoidably leads to instant 
accumulation of induced charges on the plate sides and vanishing of both an electric 
field and a current inside a metal. However, the Lifshitz formula is derived for neutral 



dp{t) 
dt 



a Q E tot (t) 



E tot (t) = E-E p (t), 



(13) 




(14) 
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plates without any nonzero surface charge densities. At the same time, the Drude model 
admits the presence of a nonzero induced current. Because of this, it is not surprising 
that the Lifshitz theory in combination with the Drude model violates the Nernst heat 
theorem for perfect crystal lattices [161 E] and was found to be in contradiction with 
several experiments [23[ EH [25], [26]. If metal plates were really infinite (as is formally 
suggested in the derivation of the Lifshitz formula) the Drude model would be applicable 
including the zero-frequency term. This, however, is an unphysical case and it cannot 
be considered as a closed system where the laws of thermodynamics must be valid. 

The above discussion uses the formulation of the Lifshitz formula ([!)) in terms of the 
imaginary frequency axis. However, direct computations using the formulation in terms 
of real frequencies show [H] that the region of sufficiently low real frequencies results in 
precisely the same contribution to the Casimir free energy as does the zero Matsubara 
frequency. As a result, all above conclusions are equally applicable to the contribution 
into the free energy from the zero-frequency term in the Matsubara formulation and to 
the equivalent contribution from low real frequencies in the formalism of real frequency 
axis. 

By contrast with the Drude model, the plasma dielectric function ( TTTj) does not lead 
to a real current of conduction electrons and does not result in accumulation of charges 
on the side surfaces of finite metal plates. The free electron plasma model in combination 
with the Lifshitz formula satisfies the Nernst heat theorem [THJ [17] . However, as was 
mentioned in the Introduction, it is in disagreement with short separation experiments 
on the measurement of the Casimir force. Below we demonstrate that the generalized 
plasma-like permittivity [36], which is in agreement with all experiments performed up 
to date also satisfies the requirements of thermodynamics. Thus, it is becoming the best 
known candidate for the adequate description of metals in the framework of the Lifshitz 
theory. 



3. Thermodynamic test for the generalized plasma-like dielectric 
permittivity 

The generalized plasma-like dielectric permittivity can be presented in the form [36] 

e(w) = l-^+A(w), (15) 

where the additional term A(u) takes into account the interband transitions of core 
electrons. Explicitly it is given by 

M") = Y,-2 y — : > (16) 

where ujj ^ are the resonant frequencies of oscillators describing the core electrons, 
gj are the respective relaxation parameters, fj are the oscillator strengths and K is the 
number of oscillators. The values of oscillator parameters for different materials can 
be found in [37]. Recently the precise determination of these parameters for Au was 
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performed in [22] • Note that the generalized plasma model does not include relaxation 
of the free conduction electrons. The latter are described by an oscillator with zero 
resonant frequency, ujq = 0, which is not contained in ( fTBT) but is explicitly included 
in (|T5j) with g Q = and f Q = u 2 . Thus, similar to the usual plasma model (fTTj) . the 
generalized plasma-like permittivity ( TL5l) . ({TBI admits only the displacement current and 
does not allow for the accumulation of charges on the sides of finite plates. Because 
of this, the generalized plasma-like permittivity is compatible with the Lifshitz formula 
which is derived for neutral plates with zero charge distributions. 

In [36] it was shown that the permittivity (|T5j) . (|T6|) precisely satisfies the Kramers- 
Kronig relations. Here we prove that the Lifshitz formula combined with the generalized 
plasma-like permittivity is in agreement with the Nernst heat theorem and thus 
withstands the thermodynamic test. 

To find the asymptotic behavior of the Casimir free energy and entropy at low 
temperature, we first present equations ([I]), ([2]) and (jT5l) . ({TBI) in terms of the following 
dimensionless parameters 

w« i . ii 
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c/(2a) is the so-called characteristic frequency of the Casimir effect. In 



where u c 

terms of new variables, the Lifshitz formula (pQ) takes the form 
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32n 2 a 3 
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{in [l -r^ u (C h y)e 
The reflection coefficients ([2]) are given by 

rrti(Ci,v) 
rTE(Ci,y) 



+ In 



ydy 
l)(2/ 2 -C 2 ) 



(18) 



(e, + l)y 2 + (e, - l)Cf + 2e iy hi{y) coth [£hi(y) 
(ei ~ l)Cf 



2j/2 + ( e , - l)Cf + 2yh l (y) coth ^(y) 



where 



Mv) = y 2 + (ei-i)Cf 



1/2 



(19) 
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The generalized plasma-like dielectric permittivity along the imaginary frequency axis 
can be presented as 

1 



+ A, 



(21) 



where 



A' 



4 = = E j 



(22) 



Using the Abel-Plana formula [5J [3] 
1 



/ \ poo poo 

E(i-^)F(i) = / F( ()d( + i / o * 



F(it) - F(-it) 



- 1 



(23) 



where -F(z) is an analytic function in the right half of the complex plane, we can 
rearrange (fTB"]) to the form 

T{a,T) = E(a)+AJ 7 (a,t). 

Here, the energy of the Casimir interaction at zero temperature is given by 

Tic 



E a) 



poo 

d( f((,y)dy 

-'C 



-v 



32vr 2 a 3 

and the function y) is defined as 

f(C,y) = [l -r 2 ?M ((,y)e~ y ] +y\n 
The thermal correction to the Casimir energy is expressed as follows: 

A^(a,T) 



r TE(C,2/)e 



ihcr [°° dt F (S Tt ) ~ F(-irt) 



32vr 2 a 3 Jo 



?2nt 



(24) 
(25) 
(26) 
(27) 
(28) 



where 

poo 

F(x) = / dyf(x,y). 

J X 

The behavior of the thermal correction (J27J), (128]) at low temperature will be the subject 
of our further consideration. 

Perturbation expansion can be performed in analogy to papers [T6] (TTJ H5] |4~6] . 
At first, we expand the reflection coefficients f|T9l) with ^ replaced by £ in powers of 
parameter a defined in ( TTTj) preserving all powers up to the fourth inclusive. The 
parameter a can be identically presented as a = A p /(47ra), where X p is the plasma 
wavelength. This means that a C 1 at all separation distances between the plates 
larger than X p . As is seen from (126]) . it is more convenient to expand the logarithmic 
functions containing the reflection coefficients in ( |26l) multiplied by the variable y. The 
results are: 

r - 1 f - \ K 2 8e y ( A 

y\n[l- r 2 TU (C,y)e »j = yln (l - e^j + a^—j - tt 2 y(eg _ 1)2 

2C 2 {2C 4 (3e« + l) 2 + 3 (e» - l) 2 y 2 [y 2 - 2C 2 - ( 2 A( 



+ a 1 



— a 



3y 2 (e» - l) 3 

8e\ 4 {2C 4 (e» + l) 2 + (e» - l) 2 y 2 [y 2 - 2( 2 - C 2 ^(C)]} 



yln l-r^ E (C,y)e" 



(ey 



dn(l 



(29) 



+ a 



e» - 1 



a: 



8y 3 e 



3 P y 



(ey - 1)' 



2y 2 -3 (e» - l) 2 C 2 -4(C) + 2/ 2 (15e 2y + 18e* - 1) 



3(e« 



8yV [- (e y - I) 2 C 2 A(Q + y 2 (e 2y + 6e y + 1) 
(& - l) 4 : 
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It is significant that these expansions do not depend on d (the thickness of the 
plates) contained in ffl"9]) . This is because the factor in the denominator of ffl"9l) . 



coth 



d 
2a 



hi(y) 



coth 




(30) 



1 — exp 
behaves asymptotically as 

d 



1 + exp + + 

1 + a V + a 2 Mi 



_d_ 

aa 



2 exp 



aa 



(31) 



when a goes to zero. Thus, this factor could only contribute exponentially small terms 
in the expansion fl29l) providing the plate thickness d is much larger than the penetration 
depth of electromagnetic oscillations into the metal [recall that 2aa = X p /(2tt)]. Under 
this condition the perturbation expansions (129]) are common for two semispaces and 
for two plates of finite thickness. We note also that terms in (1291) of order a , a and 
a 2 do not contain contributions from the core electrons. They are the same as for the 
usual free electron plasma model (fTTj) . The contributions from the core electrons are 
contained only in the terms of order a 3 and a 4 in ( T29l) . 

The parameter r defined in ( ITT)) can be identically represented as 



T 



T he 
2tt— , k B T cS = — . (32) 

Here T e g is the so-called effective temperature. For example, at a separation of a = 1 /j,m 
T e g pa 1145 K. Below we will consider the limiting case of low temperatures T <C T e g. 
The contribution from the terms of order a , a and a 2 in (j29|) into the thermal correction 
( 126|) was found in [35j H6] where the usual free electron plasma model was considered. 
This contribution is given by 

fC(3),. 1 
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where ('(z) is the Riemann zeta function. As was shown in [16], the terms in ( 1331) of 
order a and a do not contain corrections of order r n with n > 5. They contain only 
the exponentially small corrections of order exp(— 2tt/t). 

Now we deal with the terms of order a 3 and a 4 in ( |29l) which contain the 
contributions from the core electrons. From fl26|) . fl28|) and fl29l) the respective functions 
F^ 3 ^(x) and F^'(x) are given by 

y 2 dy 
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2 x e [°° (3e^ + l) 2 
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F (4) (x) =8a 4 \A(x)x 2 
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2x s 
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Calculating all integrals in (l34"j) as asymptotic expansions at small x (see Appendix 
for details) we arrive at 

A 

F®(iTt) - F^(-irt) = -2ia 3 



F (4) (irt) -F (4) (-irt) = 8ia 



2r 3 t 3 C(3)]TC j 5 : 

+ 7TT 4 t 4 + 2 

A 

3=1 



(35) 



The terms omitted in (1351) are of order r 5 . 

Substituting (135]) in (T2~Tj) and integrating with respect to t, we obtain the 
contribution to the thermal correction from the terms of order a 3 and a 4 : 
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(36) 



The total Casimir free energy computed using the generalized plasma-like 
permittivity can be now found from ( 1241 . ( 1331) and ( 1361) : 

.F(a, T) = E(a) + AF p (a, T) + AF g (a, T). (37) 

Taking into account (132|) . it can be represented in the form 

ft<(3) / T x 3 



jF(a, T) = £(a) 



7T 3 T 



167ra 3 VT e fj 



1 + 4a 



4503) T cfr 



1 + 8a + 6C(3)a 3 £ - 96C(3)a 4 ^ C^- 



96tt 2 C(5) / T 



C(3) 



off 



a: 



A' 



7T Z 



(38) 



Here one can see that the free energy calculated using the generalized plasma-like 
permittivity contains the correction of order (T/T e g) 4 not only in the terms of order a 
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and a (as in the usual plasma model) but also in the third and fourth order expansion 
terms in a. In the usual plasma model the terms of order a 3 and a 4 contain the thermal 
corrections only of order of (T/T e ff) 5 and higher [15]. To estimate the relative role of 
the additional terms arising due to the use of the generalized plasma-like permittivity, 
one can use the parameters of oscillator terms in (1221) for Au [26J. This results in 



0.272, 
0.109, 



a 
a 



200 nm, 
500 nm. 



Cj = 6.3175, '£C j 8 j = 

3=1 i=i 

From (138p it is easy to find the asymptotic behavior of the Casimir entropy 

dF{a, T) 



S(a,T) 



dT 



(39) 



(40) 



at low temperature. The result is 

3((3)k B 



S(a,T) 



S7TGr 



T 



^- <l+4a 



4tt 3 T 
13503) 7k 

160vr 2 C(5) 
C(3) 




A' 



A' 



1 + 



a 



3=1 



8a + 6C(3)a 3 ]T C^- - 96C(3)a 4 £ C j 8 j 

3=1 
K 



47r^ 




(41) 



As is seen from 

S(a,T)^0 when T -> 0, (42) 

i.e., the entropy goes to zero (and remains positive) when the temperature vanishes. 
This means that the Nernst heat theorem is satisfied and the Lifshitz theory combined 
with the generalized plasma-like dielectric permittivity withstands the thermodynamic 
test. 



4. Conclusions and discussion 



In the foregoing we have continued the elaboration of new theoretical approach to the 
thermal Casimir force based on the Lifshitz formula combined with the generalized 
plasma-like dielectric permittivity [36]. In the first part of the paper (Section 2) the 
physical reasons were presented why the Drude dielectric function is not applicable in 
the case of finite metallic plates. It was shown that for the validity of the Drude model 
the nonzero current of conduction electrons must exist, whereas the surface charge 
densities must be equal to zero. Both these conditions are shown to be violated when the 
plane wave of electromagnetic oscillations of vanishing frequency falls on a finite metal 
plate. In this case the electric field and current of conduction electrons practically 
instantaneously turn into zero. This is accompanied by the accumulation of charges 
on the sides of plates. Not only the Drude dielectric function, but also the Lifshitz 
formula are not applicable to this physical situation. On the contrary, the generalized 
plasma-like permittivity leads to only a displacement current and does not result in the 
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accumulation of surface charges. The obtained results furnish insights into the long- 
debated problem why the Lifshitz theory combined with the Drude dielectric function 
results in contradictions with thermodynamics and experiment. It becomes clear also 
why the generalized plasma-like permittivity, which does not include the relaxation 
processes of conduction electrons, is consistent with all available measurements of the 
Casimir force at both short and large separations. 

Recent paper [IT] also argues that the finite size effects of the conductors may play 
an important role in the problem of the thermal Casimir effect. This was illustrated 
in the simplified case of two wires of finite length described by the Drude model and 
interacting through the inductive coupling between Johnson currents. If the capacitive 
effects associated with the end points of the wires are not taken into account, the thermal 
interaction between the wires leads to the violation of the Nernst heat theorem. If the 
capacitive effects were taken into consideration, the agreement with thermodynamics is 
restored [Tfj . 

To conclusively establish the applicability of the generalized plasma-like 
permittivity in the theory of the thermal Casimir force between metals, in Section 3 
we have performed the thermodynamic test of this model. We have analytically found 
the asymptotic behavior of both the Casimir free energy and Casimir entropy at low 
temperature. This was done using the perturbation theory in two small parameters. 
The obtained new expressions generalize the previously known ones (found for the 
usual free electron plasma model which does not take dissipation into account). When 
the oscillator parameters describing the core electrons go to zero, the newly obtained 
expressions for the Casimir free energy and entropy go into the ones found for the usual 
plasma model. The Casimir entropy at low temperature derived using the generalized 
plasma-like permittivity is positive and takes zero value at zero temperature. Thus, the 
Nernst heat theorem is satisfied. 

To conclude, the generalized plasma-like permettivity provides a good basis in 
agreement with thermodynamics and experiment for the description of the thermal 
Casimir force between metallic plates of finite size using the standard Lifshitz theory. 
A more fundamental approach to the resolution of this problem would require, in 
accordance with Parsegian's insight [37] . the consideration from the very beginning 
of finite plates and charging of their boundary surfaces. This, however, goes far beyond 
the scope of the Lifshitz theory. 
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Appendix 



Here we derive equations ( 1351) . where functions F( 3 \x) and F^(x) are defined in (Si 

The first integral in the definition of F^(x) converges when x — > and can be 
calculated as 



A (3) (x) = /°° |^ = 2Li 3 (e-) + 2xLi 2 ( e -) - x 2 In (l 



— e 



2 S 4 

2(( 3 ) - — + — - — + 0(x 6 ) 
sv ; 2 6 48 V ; 



(a.i; 



where Li n (z) is the polylogarithm function [18]. Expanding the function A[x) defined 
in (|22|) in powers of x and using (lA.lj) we arrive at 

„4 K K 



[A(x)-1] x 2 I^\x) 



-2C(3)x 2 + %- + 2C(3) £ C 3 x 2 - 2C(3) E 

1 3=1 3=1 



A" 



A' 



A 



x 4 + 0(x E 



- E Q - 2C(3) E + 2 C(3) E <?i7i 

. Z 3=1 3=1 3=1 

From (]A.2j) . only the term proportional to x 3 contributes to (135|) . 

The second integral in the definition of F^ 3 \x) in ( 1341) also converges when x 
It can be found in the form 

y 4 (15e 2y + 18e y - 1) dy 



(A.2) 



0. 



4 3) (x) 



+ 



( e y - 1) 



15 E 



n 



71=1 



y 4 e~ ny dy 



y 4 e 2y dy 



1 4 

-16y + y + 0(x 5 ). 



/a, (1 _ e -i/)- 

This does not contribute to (1351) in the perturbation orders under consideration. 
The third integral in F^ 3 \x), 

e~ y dy 



(A.3) 



4 3) (*) 



x _ e _ y -ln(l-c-), (A.4) 

diverges when x goes to zero. It should be, however, multiplied by [A(x) + 2] x 4 with 
the result 



[A(x)+2] x 4 4 3) (x) 



A 



A 



EQ + 2 x 4 In x + J2 C i S J x5 In x + O (x 5 ) . (A.5) 

V3=l / 3=1 

Only the first term on the right-hand side of (1A.5I) contributes to (1351) . This contribution 
is simply found when taken into account that 

ln(iz) - ln(-iz) = ivr. (A.6) 

The last, fourth, integral in the definition of F^(x), also diverges when x goes to 
zero. It can be identically represented in the form 



4 3) (*) 



e- y (3 + e~ y y 



poo p y p 
dy = 16 / — dy + 16 / 

Jx V a Jx 



x y 2 (l- e~ y ) 
00 P -y \g poo e -y 



00 p-y 



roc e -y iy roc e -y 4i r, 



x y 

41 r 00 e~ y 



x r 



-dy 



y 



dy 



+ 



00 p-y 



x y* 



(3 + e^) 2 16 16 9 19 Aly 



'I- e-y) 



6 



60 



dy. 



(A.7) 
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The last integral on the right-hand side of (1A.7j) converges when x — > 0. It does not 
contribute to the perturbation orders of our interest after the multiplication by x 5 . As 
a result, we arrive at 



x G lf\x) 



x 



16IX-4, x) + 16r(-3, x) + 9r(-2, x) 



19. 



41. 



+-r(-i,x) + -r(-o,x) 



+ 0(x 6 



where T(n,x) is the incomplete gamma function |49j. Using the identity [19 

n-1 

r(0,x) -e~ x £(-ir 



-n, x) 



tv. 



ml 



m=0 



X 



m+1 



where n — 1, 2, ... , and the asymptotic relation 



x 



X 



T(0, x) = -7 - lnx + x - — + — + 0(x 4 ) 

4 lo 

with Euler's constant 7 = 0.577216, we finally obtain 



x 



4x z + — + 0(x 5 ) 



(Ai 



(A.9) 



(A.10) 



(A.11J 



This evidently does not contribute to (1351) . 

As a result, by using (1A.2j) and (1A.5[) we obtain the first equation in (1351) . 

Now we consider the derivation of the second equation in (1351) containing the 
function F^(x) defined in fl34l) . The first integral in 



4 4) (x) 



y 3 e y 
'1 - e-y) 



converges when x goes to zero. It is calculated in analogy to 
following expansion is obtained: 

K K 

A(x)x 2 l[ A \x) = 6C(3) C 3 x 2 - 6C(3) ^ C 3 5 3 x 3 



(A.12) 

31) . As a result the 



1 



A" 



A 



A 



x 4 + 0(x 5 



(A.13) 



EC J -6C(3)E^ + 6C(3)E^- 
j=i j=i j=i 

The second term on the right-hand side of (1A.13I) contributes to (|35l) . 

The second integral in F^ 4 ^(x) also converges when x goes to zero. It can be 
calculated similar to li 3 \x) in (1A.2|) and does not contain odd powers of x lower than 
x 5 . Thus, this integral does not contribute to ( |35i) . 

The third and fifth integrals in the definition of F^ 4 ^(x), (!34j) . diverge in the limit 
x — > 0. However, by calculating them similar to I^\x) in (1A.7h f fA.10f) and multiplying 
the results by x 6 and x 8 , respectively, we find that both these integrals do not contribute 
to d35}. 

The fourth integral in F^(x) can be calculated as follows: 

coo ye~ y dy 



r(4), 



x = 



(l - e-yy 



In 1 



+ 



xe 



1 



(A.14) 
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It diverges when x goes to zero. After multiplication by x 4 one obtains 



x 4 I^\x) 



X 



x 4 In a; + 0(x 6 



(A.15) 



The second term on the right-hand side of ( 1A.15I) contributes to ( 1351) . 
Using (1A.13j) and (1A.15j) we obtain the second equation in (|35|) . 
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